y- vibrational states in superheavy nuclei 
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Recent experimental advances have made it possible to study excited structure in superheavy nuclei. The 
observed states have often been interpreted as quasi-particle excitations. We show that in superheavy nuclei 
collective vibrations systematically appear as low-energy excitation modes. By using the microscopic Triaxial 
Projected Shell Model, we make a detailed prediction on /-vibrational states and their El transition probabilities 
to the ground state band in Fermium and Nobelium isotopes where active structure research is going on, and in 
270 Ds, the heaviest isotope where decay data have been obtained for the ground-state and for an isomeric state. 



PACS numbers: 21.10.Re, 21.60.Cs, 23.20.Lv, 27.90.+b 

One of the important predictions in nuclear physics is the 
emergence of a region of long-lived superheavy elements be- 
yond the actinides, the so-called 'island of stability'. The 
question concerns the precise location of the next closed nu- 
cleon shells beyond Z = 82 and N = 126. To reach the is- 
land, much of the experimental effort has been focused on 
the direct creation of superheavy elements. In recent years, 
progress has also been made in structure studies for nuclei 
beyond Fermium, thanks to the development of detector sys- 
tems for decay and in-beam studies using recoil separators and 
heavy ion fusion reactions |Ul 01- It has been suggested that 
by studying the transfermium nuclei, in particular their ex- 
cited structure, one can gain useful information on relevant 
single-particle states [3], which is key to locating the island. 

The nuclei of our interest, the Fm (Z = 100) and No (Z = 
102) isotopes, belong to the heaviest mass region where struc- 
ture can currently be studied experimentally. The yrast prop- 
erty of these nuclei shows that they are generally good ro- 
tors. Rotational behavior of some of these yrast bands has 
been successfully reproduced by several models (see, for ex- 
ample, Refs. J3, El IS 01)- The discussion on excited con- 
figurations so far have been focused on quasi-particle exci- 
tations ISJ |Sl LLOl . On the other hand, a deformed rotor 
can, according to the collective model, undergo dynamical os- 
cillations around the equilibrium shape, resulting in various 
low-lying collective vibrational states. Ellipsoidal oscillation 
of the shape is well known as / vibration B12I1 . It is thus nat- 
ural to consider / vibrational states in superheavy nuclei, and 
in fact, this excitation mode has begun to draw one's atten- 
tion lfl3l [l4ll . Knowledge on vibrational states in superheavy 
nuclei is particularly useful for this less known mass region 
because of the interpretation of the observed low-lying spec- 
troscopy. 

Early calculation of /-vibrational states in heavy nuclei 
was performed by Marshalek and Rasmussen IU5I1 using the 
qausi-boson approximation, and by Bes et al. lfl6Tl using 
the quadrupole-plus-pairing model based on deformed Nils- 
son states 11711 . Modern treatment of y vibration includes 



the Tamm-Dancoff method, the random phase approximation 
iH, and others EES. In Ref. Ell, a shell-model-type 
method for describing /-vibrational states was introduced, 
which is based on the Triaxial Projected Shell Model B22I1 . an 
generalized version of the original Projected Shell Model 02311 
by extending it to a triaxially deformed basis. It was shown 
112 ill that by performing diagonalization in a basis constructed 
with exact three-dimensional angular-momentum-projection 
on triaxially deformed states, it is feasible to describe /- 
vibrational states in a shell model framework. In this way, one 
can achieve a unified treatment of ground-state band (g band) 
and multiphonon /-vibrational bands (/band) in one calcula- 
tion, and the results can be quantitatively compared with data. 
The underlying physical picture of generating /-vibration in 
deformed nuclei is analogous to the classical picture of Davy- 
dov and Filippov [24]. Subsequent papers 1251 l26l 12711 studied 
electromagnetic transitions by using wave functions obtained 
from the shell model diagonalization. 

The Projected Shell Model B23I1 is a shell model that uses 
deformed bases and the projection technique. In the present 
calculation, we use the triaxially-deformed Nilsson plus BCS 
basis |<t>). The Nilsson potential is 
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where Ho is the spherical single-particle Hamiltonian with 
inclusion of appropriate spin-orbit forces parameterized by 
Bengtsson and Ragnarsson B28I1 . The axial and triaxial parts 
of the Nilsson potential in Eq. ([TJ contain the deformation 
parameters e and e' respectively, which are related to the con- 
ventional triaxiality parameter by / = tan _1 (e'/e). Pairing 
correlation in the Nilsson states is taken into account by a 
standard BCS calculation, within a model space of three ma- 
jor shells for each kind of nucleon (N = 5, 6, 7 for neutrons 
and N = 4,5,6 for protons). 

The rotational invariant two-body Hamiltonian 
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FIG. 1 : Calculated energies of the g band (solid lines) and y band 
(dashed lines) in 256 Fm as functions of triaxiality parameter e' for 
various angular momenta. Each curve is labeled by Ik- The experi- 
mental g band (filled triangle-up) and possible 7 band (filled triangle- 
down) are best reproduced with e' = 0. 105. Data are taken from Ref. 



is diagonalized in the three-dimensional angular-momentum- 
projected basis {P' MK |<I>) ,0 < K < /}, where P' MK is the 
angular-momentum-projector lfl9ll . The wave function thus 
takes the form 



0<K<I 



(3) 



where a specifies the states with the same angular momen- 
tum I. The wave function (0) is explicitly written as a super- 
position of projected K states. The two-body forces in Eq. 
(O are quadrupole-quadrupole (£20, monopole-pairing, and 
quadrupole-pairing interaction, respectively. The strengths of 
the monopole and quadrupole pairing forces are given respec- 
tively by Gm and Gq in Eq. (|2]), where 
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with " + " for protons and " — " for neutrons. We use 
G\ = 21.24, G2 = 13.86, and assume the quadrupole-pairing 
strength in Eq. (O to be Gq = 0.1 3Gm, which are found to 
be appropriate for this mass region J9l |29(]. The QQ-foice 
strength % 1S determined such that it holds a self-consistent 
relation with the quadrupole deformation e [23]. 

The nucleus 256 Fm is taken as our discussion example be- 
cause it is perhaps the only nucleus in this mass region where 
some observed states were thought to be y vibrational states 
i33tl . Figure 1 shows the calculated energies as functions 
of triaxiality parameter e' for various angular momenta. Al- 
though this figure looks similar to the one shown in the semi- 
nal paper of Davydov and Filippov |24J], it is obtained from a 



TABLE I: Axial and triaxial deformation parameters used in the cal- 
culation. 
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fully microscopic theory. Unlike the asymmetric rotor model, 
our spectrum depends not only on the deformation parameters 
but also microscopically on the detailed shell filling. In Fig. 1, 
each curve is labeled by I a . It turns out (see later discussions) 
that for the nuclei studied in the present paper, K mixing is 
very small. Therefore, a ~ K, and we can practically use Ik 
to label the states. One sees that, for the g band (K = 0) in 
256 Fm, the energies are nearly flat as e' varies, and the val- 
ues remain very close to those at zero triaxiality. Thus, we 
can conclude that the triaxial basis has no significant effect on 
the g band in such a superheavy rotor, and one can practically 
describe g band using an axially deformed basis. 

On the other hand, it has a large effect on excited bands 
with K (second and higher excited bands are not shown 
in this figure, but will be discussed in Fig. 3). Their excita- 
tion energies are indeed very high at e' « 0, but come down 
quickly as the triaxiality in the basis increases. At e' = 0.105, 
the first excited band {K = 2) nicely reproduces the observed 
2 + band in 256 Fm H33h . In our calculation, the triaxial param- 
eter e' serves as a free parameter adjusted to reproduce the 2 + 
bandhead. It should be noted that the excited bands studied in 
this paper are obtained by introducing triaxiality in the basis 
(quasiparticle vacuum). They are collective excitations, not 
quasi-particle excitations. 

The above results can be further understood by studying 
the K mixing coefficients for each projected K state in the 
total wave function of Eq. (0. It is found that for 256 Fm (and 
for all nuclei studied in this paper), K mixing is negligibly 
small. States in the g band are essentially the projected K — 
state for any e'. The basis triaxiality does not influence the 
g band result, and the rotational behavior of the g band can 
be understood by using a simple axialry-deformed rotor. The 
excited bands are also built by rather pure projected K states. 
For example, the first excited band with the bandhead spin 
/ = 2 is mainly the projected K = 2 state (labeled as Ii). The 
fact of very small K mixing in these nuclei sets up a favorable 
condition for the occurrence of K isomeric states sea. 

In Fig. 2, we present results for even-even Fm and No iso- 
topes with neutron number from N = 150 to 156. This is the 
mass region where active structure study is currently going on. 
In Fig. 2, theoretical results for both g and 7 bands are pre- 
dicted up to / = 10. The axial and triaxial deformation param- 
eters used in the present calculations are listed in Table I. We 
note that deformation in nuclei is a model-dependent concept. 
Our deformations are input parameters for the deformed basis, 
and in principle, it is not required that the numbers in Table 
I are exactly the same as nuclear deformations suggested by 
other models. Nevertheless, it turns out that our employed ax- 
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FIG. 2: Calculated energies for g bands (open circles) and / bands (open squares). Available experimental data (filled triangle-up for g bands 
and triangle-down for 7 bands) are also shown for comparison. Data are taken from Refs. 03 ( 250 Fm), 01 ( 252 Fm), 01 ( 254 Fm), 01 
( 256 Fm), 01 ( 252 No), and Offl ( 254 No). 



ial deformation parameters are very close to those calculated 
in Refs. USUI, and follow the same variation trend along 
an isotopic chain as predicted by other models (for example, 
the most deformed isotope has the neutron number 152 and 
a decreasing trend for heavier isotopes is expected). A strict 
test for using these parameters will be whether the calculation 
can describe all observables, the most relevant one being the 
B(E2) value (see discussions below). The triaxial parameter 
e' = 0. 105 that gives the correct position of the experimentally 
observed excited band in 256 Fm 13311 corresponds to 7 = 25° 
in terms of the usual gamma parameter, if one uses as a rough 
estimate 7= tan -1 (e'/e). Note that there are no existing data 
that are firmly assigned to be 7-vibrational states. The excited 
band starting from I n = 2 + in 256 Fm IB3I1 was assumed by 
Hall et al. to be 7-vibrational states. Now our results strongly 
support the interpretation of Hall's 2 + band as 7-vibrational 
band. For those nuclei that have no 7 band data to compare 
with, we simply employ e' = 0.1 in the calculation. Over- 
all, we predict 7-vibrational bands in these nuclei, with the 
rotational behavior very similar to that of g bands. With the 
present set of parameters, we find that the Ik = 1i bandhead 
energy lies low, generally at 0.6 — 0.8 MeV above the ground 
state. 



The type of a-decay experiment [ 39] makes the identifi- 
cation of excited states in the heaviest nuclei possible. We 
take an example from the heaviest isotopes and perform cal- 
culations for 270 Ds. With the same interaction strengths as 
used for the Fm and No isotopes and deformation parameters 
e = 0.225 and e! — 0. 1, we obtain the sequence of ground state 
rotational band energy E(I K = 2 ) = 0.048, £(4 ) = 0.159, 
£(6 ) = 0.332, £(8 ) = 0.568, and £(10 ) = 0.863 (all in 
MeV), in a good agreement with the experimental estimation 
lt39Tl . Furthermore, we predict very low-lying collective 7 vi- 
brations, with the bandhead energy £(22) = 0.565 MeV. The 
results will be presented later in the discussion on 7 phonon 
states. We note that at present, there is no better way to deter- 



mine the triaxiality for 270 Ds, and using e' = 0. 1 for the whole 
mass region is a natural choice. What we want to emphasize 
is that with this assumption, there is a low-lying 7 band in this 
truly superheavy nucleus. 

Recently, a great deal of experimental effort has been made 
to understand low-lying excited structure in transfermium nu- 
clei. The discussion so far has been focused on quasi-particle 
excitations only (§, 0, Qji !!]■ Our present results suggest 
that there must exist an important part of low-energy collec- 
tive vibrations, in particular at the energy region below 1 MeV. 
It is easy to make the following estimation. The pairing en- 
ergy gaps from our BCS calculation for 256 Fm are A w 0.5 
MeV for neutrons and A sa 0.8 MeV for protons, which are 
needed amount to correctly reproduce the rotational behavior 
(e.g. moment of inertia). Thus, energy of a 2-quasiparticle 
state should be greater than 2A, i.e. 1 MeV. Therefore, we 
expect that a 2-quasiparticle excitation energy is generally 
higher, which lies at or above 1 MeV in this mass region. 
This is a useful criterion to identify a low-lying excited state 
as a collective state rather than a quasiparticle excitation. Of 
course, certain residual interactions may push the quasiparti- 
cle configurations down. Very recently, possible occurrence 
of low-lying alternative parity bands in superheavy nuclei due 
to octupole correlations has been suggested [i40ll . 

A detailed prediction for 7 band depends on pairing 
strengths which are parameters in the model. The parame- 
ters G\ and G2 in (0|i are adjusted in accordance with the size 
of single-particle space, and are chosen so that they can give 
correct rotational sequences for the superheavy mass region 
(as presented in Fig. 2). Figure 3 shows the deviation from 
the current prediction if the pairing strengths Gm, and thus 
Gq (since Gq is proportional to Gm), are allowed to vary by 
±10%. It is observed that the curvature of the curves, i.e. the 
rotational frequency co = AE/AI, increases with increasing 
pairing strength. Pairing also shifts the 7 bandhead; a stronger 
pairing leads to a higher bandhead. In Fig. 3, about 20% of 



4 



q3 
c 
LU 



1.5 



0.5 



I 1 I 1 

Gyx0.9 


' I ' I 
/ 

/ 

y 


G M x1.1 


y y 
y / 
y ^y 




^ y 
y , 
y / 
y /..-'*' 
y /\.. 


| , L 


I , I 



4 6 
Spin 



10 



FIG. 3: Comparison of calculated g band (lower group of curves) and 
7 band (upper group of curves) for three sets of pairing parameters: 
the original pairing (solid curves, which reproduce the data, see Fig. 
2), the one multiplied by a 0.9 factor (dotted curves), and the one 
multiplied by a 1.1 factor. 256 Fm is taken as the example. 
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FIG. 4: Predicted multi-phonon y-vibrational spectra in 256 Fm and 
270 Ds. Available experimental data of 256 Fm 1 33] for the g band and 
7 band are also shown. 



deviation is seen for the y bandhead when pairing changes by 
10%. However, even with this amount of uncertainty in pair- 
ing, the predicted y bandhead energy is still within 1 MeV of 
excitation, and thus the conclusion for the occurrence of low- 
energy 7 vibrations remains valid. 

Multi-phonon y bands are rotational bands built on top of 
7 vibration classified by phonons. In Fig. 4, we plot all the 
states for spins 7 < 10 obtained by diagonalization for 256 Fm 
and 270 Ds. One sees that the excited states are clearly grouped 

according to K — 0,2,4, We identify the first excited band 

(K = 2) as the y band, the second excited band (K = 4) as 
the 27 phonon band, the third excited band (K — 6) as the 3y 
phonon band, and so on, with presence of strong anharmonic- 
ity in vibration. We have calculated all eight nuclei discussed 
in Fig. 2, and obtained similar patterns as for 256 Fm. In partic- 
ular, we predict 27 phonon bands for them with the bandhead 
energy at about 1.8 MeV. 27 phonon states were observed in 
the rare earth nuclei 166 Er and 168 Er. To compare with those 
27 phonon bandhead energies, 2.029 MeV in 166 Er 14J and 
2.056 MeV in 168 Er [42], 27 phonon states in superheavy nu- 
clei lie lower. No 3y phonon states have yet been seen ex- 
perimentally in any known examples of nuclear spectroscopy. 
According to our calculations, they should appear at about 3 



MeV in Fm and No isotopes. To compare with 256 Fm, a more 
compressed multi-phonon 7 vibrational spectrum is seen for 
270 Ds, again with strong anharmonicity. 

The wave functions obtained after diagonalization of the 
Hamiltonian are used to calculate the electric quadrupole tran- 
sition probabilities 



B(E2:(I h Ki)^(If,K f )) 



27,; +1 



between an initial state (Ii,K() and a final states (If,Kf). The 
explicit expression for the reduced matrix element in the pro- 
jected basis can be found in Ref. Il25ll . Note that we now use 
K instead of a to specify states with the same angular mo- 
mentum I to keep the familiar convention. In the calculation, 
we use the standard effective charges of 1 .5e for protons and 
0.5e for neutrons. 

In Table II, we list calculated B(E2) values within the g 
bands, the inter-band linking transitions between 7 bands and 
g bands, and between 27 bands and 7 bands. It is found that 
for each of the nuclei, the g band B(E2) values correspond 
to a rather constant transition quadrupole moment, reflecting 
the fact that these systems are good rotors. The inter-band 
transitions are on average by two to three orders of magnitude 
smaller than the in-band transitions within the g bands. These 
predicted B(E2) values may help in determining the structure 
of low-lying states. 

To summarize, we have applied the Triaxial Projected Shell 
Model to Fm and No isotopes as well as 270 Ds to predict 7 
vibrational states in superheavy nuclei. Shell model diag- 
onalization is carried out in a angular-momentum-projected 
triaxially-deformed basis. It is found that the calculation si- 
multaneously leads to a consistent description of ground state 
band and multi-phonon 7 bands in these nuclei. The physics 
of the bands is discussed in terms of K mixing, suggesting a 
microscopic connection between the excited vibrational states 
and the nuclear ground state properties in the heaviest nuclei 
where structure study can be performed experimentally. In- 
and inter-band B(E2) values are predicted. This work calls 
for attention on collective excitations in the low-energy region 
where active structure studies are currently carried out. 
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TABLE II: Calculated in-band B(E2,Iq — > (/ — 2)q) values for g bands, inter-band transitions from y bands to g bands B{E2,lj_ — * ^o). and 
from 2y bands to y bands B(E2,li t — > /j). Numbers are given in W.u.. 
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